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Abstract 

We present some techniques which have been developed recently or in the recent past to 
compute Feynman graphs beyond one-loop order. These techniques are useful to compute 
the three-loop splitting functions in QCD and to obtain the complete second order QED 
corrections to Bhabha scattering. 



1 Introduction 

Since the discovery of asymptotic freedom in non-Abelian gauge field theories, like Quantum 
Chromo Dynamics (QCD), many perturbative calculations have been performed to hadron- 
hadron and semi-leptonic processes. The most reactions are computed up to next-to-leading 
order (NLO) only, except for some semi-leptonic processes. The latter is due to the simplicity of 
the Born approximation where the basic process is given by the interaction of the intermediate 
vector boson with a quark. Higher order corrections to the other reactions are still missing be- 
cause of the very complicated Feynman and phase space integrals which arise in the calculations. 
Here we will discuss some new methods which maybe will enable us to compute the QCD and 
also QED corrections beyond the NLO level. 
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2 Asymptotic Expansions 



In this section we concentrate on the computation of the coefficient functions and anomalous 
dimensions which show up in deep inelastic lepton-hadron scattering. The relevant quantity in 
this process is the structure function defined by 

F(x,Q 2 ) = -^Jd*ze^(p\[J(z),J(0)}\p} q 2 = -Q 2 <0 x = ^-. (1) 

In the commutator above J(z) represents the electro-weak current where we have suppressed 
Lorentz indices for convenience. In the Bjorken limit i.e. Q 2 — ■> oo and x = fixed the light cone 
region dominates the integrand so that one can make an Operator Product Expansion (OPE) 



[J(z),J(0)] z CoY.C N (z 2 V 2 )0 N (^,0). 



(2) 



N 



This expansion also holds for the time ordered product corresponding to forward Compton 
scattering 



T(x,Q 2 ) 



d 4 ze^(p\[J(z),J(0)]\p), 



An 2 

ImT(x,Q 2 ) = 7tF(x,Q 2 ). (3) 
Let us assume that we can write an unsubtracted dispersion relation in the variable v = p ■ q 

t(v, Q') = J°°_ i"'—^r~r~ . (4) 



>Q 2 /2 



v — V 

2 



so that after substitution of the variables v — v' — , 



we can derive the following relation 



Jo x x — x' 



N=0 



where the operator matrix element and the coefficient function are given by 

A( N \^) = (p\O N ( f , 2 ,0)\p), 

and 

CW (91^ = J 'd 4 ze^ z C N (z 2 ii 2 ), 
respectively. Both quantities satisfy a renormalization group equation 

#)(^ 2 ) = 0, 



(5) 
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Figure 1: Reduction of a two-loop Compton graph into a self energy diagram. 



In the equations above the quantities /3(g) and j^^g) represent the beta- function and the 
anomalous dimension respectively. The latter determine the Q 2 -evolution of the structure func- 
tion F(x, Q 2 ) which can be measured in experiment and provides us with one of the tests of 
perturbative QCD. If one computes the coefficient functions in the conventional way one en- 
counters phase space integrals and loop integrals (see e.g. U). However if one wants to compute 
the coefficient function beyond order g A it is more convenient to try another method which is 
explained in 0. Taking 0g-theory as an example one can expand the propagator in T(x, Q 2 ) (H) 
given by l/(k — p) 2 as follows 



T(x,Q 2 )=T = (-ig) 2 J 



d n k i 
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(2tt)» (k 2 ) 2 (k - p) 2 (k + qf 



= -iff I**- i2k - p)N (9) 
9 tJ (2iry(k 2 ? +N (k + q) 2 ' {J> 

where we have used n-dimensional regularization to regularize the collinear (C) and ultraviolet 
(UV) singularities. The asymptotic expansion of the propagator converts the Compton amplitude 
into a self energy type of integral. The latter is much easier to compute than expressions for 
box- or triangle graphs. Another feature is that this asymptotic expansion transforms the C- 
divergence in T at n = 6 into an UV singularity appearing in the expression for the self energy 
integral. Therefore the operator matrix element in Eq. (|5|), containing the collinear divergence 
denoted by 1/ec, will be replaced by the corresponding operator renormalization constant Z n 




/ \£uv / V A* / 

where Z q n depends on the UV singularity denoted by 1/euv- This transformation leaves the 
coefficient function unaltered. A straightforward calculation yields 

2 vr"/ 2 r(n/2 - l)I>/2 - 3) - (2p ■ g\ N T(4 + N - n/2) (Q 2 ^ 

w 9 (2tt)" r(n-4) tX Q 2 ) r(3 + A0 UV ' 1 } 

One can expand the gamma-functions above around e = n — 6 and obtain 



F ^ 2p-q 



N 



647r 3 ^ n V Q 2 I Lat + i n + 2 
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Figure 2: Diagrammatic representation of Eq. (|16|). The symbol || means that the corresponding 
propagator appears twice. 



n 2 N 1 



In order g 2 Eq. (p70|) can be expressed as 



2 oo 



7 ^ f-- 7 y, -h.47r-ln^^ * 

e /r 



(13) 



so that one can read 7^ and c^ N \ If the method is extended up to order g 4 one encounters 
graphs where at least three propagators, carrying the momentum p (see the solid line in Fig. 
(0)), have to be expanded This procedure leads to triple sums which are very difficult to unravel 

( 2p-k 1 Y ( 2p-k 2 \ j ( 2p-(k 1 -k 2 -q) \ S 

To avoid this one has to remove one or more propagators, carrying the momentum p, before 
doing the expansion. First one can adopt the method 

1. Integration by Parts 

This trick has been successfully applied to compute self energy graphs ||. An example is 

d n h f d n k 2 1 

(15) 



(2tt)» J (2tt)« k\k\(k x + q)\k 2 + qfik^ - k 2 f ' 

The integral above appears in the partial integration w.r.t. the momentum k\ of the 
expression below (see Fig. (0)) 

[d^rcPk^Jl (h - k 2 )^ 



(2tt)» 7 (2tt)« 9A;i, p fc?fci(fci + q)\k 2 + g) 2 ^ - fc 2 ) 

Hence / in Eq. ([HJ) can be reduced to a set of integrals each containing one propagator 
less than in the original expression. 
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Figure 3: Difference Equations. 
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Figure 4: Building Blocks. 



2. Recursion Relations (Difference Equations) 

Another method is developed in ||. Using the method of "Difference Equations" one is 
able to knock out one propagator (see the example in Fig. (H)). One can show that all 
graphs can be reduced to two building blocks only (see Fig. @). The latter contain one 
propagator carrying the momentum p so that one only obtains a single sum which is easy 
to perform. At this moment it is not clear whether the method can be extended beyond 
order g A (two- loop). 

In higher order it is only possible to obtain the integral over the structure function in Eq. 
(H) for finite moments iV (see [|]]). This is achieved by 



d N T 



In our example Eq. flj^) this procedure leads to the result 



2^ N\ % \'^ q S N Z n 



)2W 



Q 2 



2 N N\(-igf 



dPm ■ ■ ■ dp m 
which is again of the self energy type. 



• ■ h 



(2tt) w (k 2 ) N+2 (k + qf 



(17) 
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3 Mellin Barnes Techniques 

A method which turns out to be very useful to compute three-point and four-point functions 
even up to two- loop level is given by the Mellin-Barnes technique ||. A simple Mellin-Barnes 
transformation takes the form 



1 rioo 

(A 1 + A 2 )-» = — A\A- 2 »-° 

LlXl J-ioo 



(19) 
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with \arg(Ai) — arg(A 2 )\ < 7r and the poles in a are located on the real axis. To illustrate this 
technique we take as an example the three-point function where the propagators are raised to 
an arbitrary power 



J(n; z/i, v 2 , 1/3) = J d n k 



1 



+ k) 2 - m 2 }" 1 [(q 2 + k) 2 - mjp [(q 3 + k) 2 - rri 2 ]^ 
Zi v% ~ n , 

n,r(z/ l ) 



»^ r ^ , ? /2) / n**«r 1 *a-5>i 



(20) 



where Pi = q-s — q 2 , p 2 = qi — 93 and P3 = q 2 — qi- Next we apply the double Mellin Barnes 
integral transformation given by 

1 1 r ico du f io ° dv T(a + u + v)T(—u)T(—v) 



(X + Y + Z) a Z a J^oo 2m J-ioo 2m T(a) 



1 



(21) 



so that the integral in Eq. fl20|) can be written as 



H — r — n^-n/2+u+vn-un-v) ^4 ^4 , (22 ) 

J-ioo 2m J -too 2m, KZ f 1 \aiplj \a 2 plj K J 
where we have put for simplicity m 2 = 0. Using 

tt (\ X n ^ a b c r(i + a)r(i + 6)r(i + c) , O0 > 

/ daidil - > a,' a?a,a§ = ' ' , (23) 



the three-point function reads 
J(rc; ux, u 2 , u 3 ) = 



7r"/ 2 r 



i_„ T 00 dw r°° r(n/2 — - z/ 3 )r(n/2 — ^ 2 - ^3)r(i/3 + u + f) 
ioo 27tz J-ioo 2m Hi T(ui)T(n — J2j v j) 



T^k + u + v- n/2)Y{-u)Y{-v) (^) " " (p 3 ) n/2 ~^ • 



(24) 



A special case is z/i + z/ 2 + z/ 3 = n (Uniqueness Condition). After application of the residue 
theorem to expression (|4]) one obtains 

J(n; Vl , V2 , „) = f[ r( "ff7 J> (p|)"'-" /2 . (25) 

j=i 1 
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Figure 5: Two-loop triangle graphs. 




Figure 6: Two-loop box graphs. 



The Mellin Barnes techniques has been successfully applied to compute two-loop graphs. Ex- 
amples are the three-point functions in Fig. |5] and the four point functions in Fig. ^| which are 
computed in || for p\ ^ 0. Recent progress has been made for the two- loop box graphs in Fig. |6| 
which has been computed for all external momenta on-shell (pf = 0) ||, 0. This is very useful 
for applications to the second order corrections in Bhabha scattering and in di-jet production. 



4 Negative Dimension Approach 

This method has been recently developed in || || to express one-loop integrals containing 
different masses into special functions. To illustrate the technique we take the two-point function 
as an example. The latter is given by the integral 

r d n k 1 

r 2 {^ 2) q\MlMl) = j inn/2A „ 1A „ 2 , (26) 

Using the Schwinger representation for the propagator 

1 (-1)" 

A? rfo) 

expression ( ^6|) becomes equal to 



dxiX^ exp(xiAi) , (27) 



% = II WtV /„ *r 1 / -^r 2 exp I £ XjAj I . (28) 



\ T(z/j) Jo 1 1 J m n l 2 



The expression above can be evaluated in two different ways. The first way proceeds by expanding 
the exponents in a power series. This yields 



2 -r (-I)" 1 r , f d n k ~ (x 1 A 1 ) n - ~ (x 2 A 2 ] 



r > - n r( ^y y^s, ^ £ n 2! 



"2 
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2 (—D^i r°° 00 00 2 r nj 

II / dxi E E ~n2, q\ Ml M\) J] • (29) 



J r(^) 7o 

Using the identity 

(_!)*. 



T(vi) Jo 

we obtain 



cfa^ 1 " 1- ™' =6„ i+n:jt o, (30) 



JJ = JJfa, u 2 , q\ Ml Ml) J] * r . (31) 

i=l 1 1 1 



There is a second way to compute Eq. (p8[). First one shifts the momentum which becomes 
k = k' — q x 2 /(xi + x 2 ). Subsequently one uses the following identities 

d n k' ,„si .. „ r d n k 



lit 



n/2 



1 r ri n h 1 

(fc' 2 ) =/M tww /^K) = ^. (32) 

Notice that since Z > one obtains n < 0. This the reason for the name "Negative Dimension 
Approach" . The expression in Eq. ( p8|) can now be written as 

= J f ^j^r^P {^-/) «* (- t Xl Mf) . (33) 

The integrand can be written as 



1 / ^1^2 2 



(xi + x 2 ) n/2 \Xi + x 2 



q 2 ^ exp ^— XjM 2 ^j 



tk^ + x2r l - n/2 (^ 2 Y E t^Mt^iMIIl . (34) 

2=0 m=0 m - 



Using the identities 



M M oo a M 

\N r/» r i n FT \ ~> ^ik^ik 



£>A.r =r(iv + i)n E ^r 1 with T,^ = n, (35) 

i=0 fc=li fe =0 %k ' k=l 



one can write 



i=i 1 K p i) Jo 1=0 L[ Pl =o P2 =aVr- 
E E ~ j Wlth Pi+P2 = -l~n/2. (36) 

mi =0m 2 =0 m l ! m2! 

Equating Eqs. ( PT| ) and (BB) one obtains the result 

r( „ „ „2 M 2 M 2, v r(i - ; - n/2)r(i - n)r(i - ^ 

" ,^,r(i + ;)r(i + Pl )r(i + p 2 )r(i + mi) r(i + m 2 ) 



/_M 2 \ mi ( '-M 2 \ m2 
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with 



Vi+V2 = -I - n/2 p 1 = l — v 1 — miP2 = l-v 2 -m2. (3£ 



Choosing two independent summation indices out of five i.e. £,pi,p 2 , mi, m 2 one can express 1% 
into the following special functions 0. given by the Hyper-geometric Functions 2F1, 3F2, the 
Appell Functions Fi (i = 1 — 4), the Horn Function H2 and the Kampe de Feriet Functions 
Si (i — 1,2). In the case of the two-point function in Eq. fl26|) there are eight possibilties to 
express the integral into the special functions above depending on the convergence of the sums 
in a specific kinematic region. If we choose \Mf/q 2 \ < 1 one has to take mi,m 2 . Using the 
following identity 

£<£_!!) = (_!,. rOzi) (39) 

T(z) r(l — z + n) 

the expression in Eq. ^ can be written as 

sew. £ £ Z^iVTV"^ 

mi =om 2 =o 1 \ n l l + l ~ v i- u 2~ mi-m 2 ) 



2 



r(l + mi)r(l + m 2 ) \ Q 2 J \ Q 

. / 1 ^/2^2^/2- i / 1 -v 2 r W 2 - "i)r(n/2 - i/ 2 )r(i^ + i/ 2 - n/2) 



r(z/i)r(z/ 2 )r(n - ^ - z/ 2 ) 

M 2 M 2 

F 4 (l + z/ x + z/ 2 - n, i/i + u 2 - n/2, 1 + v x - n/2, 1 + u 2 - n/2, -i, -|-) . (40) 



Here the fourth Appell Function is given by 

FJa 3 7 V x y) - T na + m + n)T(B + m + n) F( 7 ) F(V) 

41 ^ ,7,7, ,yj ^=0 r(«) r(/3) r( 7 + m)r( 7 ' + n) 

x m y n 



r(m+l)r(n+l) 



(41) 



We conclude that using the method above one can find different representations for an one- 
loop integral in different kinematic regions depending on the radius of convergence of the ratios 
between the scales g 2 ,M 2 ,M|. It is unclear whether this method, which seems to us very 
complicate, can be used to compute two-loop integrals. 
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